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^ ■ ABSTRACT 

O: 

. Multiwavelength study of extended astronomical objects requires combining 

images from instruments with differing point-spread functions (PSFs). We de- 
^ . scribe the construction of convolution kernels that allow one to generate (mul- 

^ . tiwavelength) images with a common PSF, thus preserving the colors of the 

^ . astronomical sources. We generate convolution kernels for the cameras of the 

Spitzer Space Telescope, Herschel Space Observatory, Galaxy Evolution Explorer 
(GALEX), Wide-field Infrared Survey Explorer (WISE), ground-based optical 
^ . telescopes (Moffat functions and sum of Gaussians), and Gaussian PSFs. Ker- 

I— I, nels for other telescopes including IRAS, AKARI, and Planck, are currently being 

^ ! constructed. These kernels allow the study of the spectral energy distribution 

1^ ! (SED) of extended objects, preserving the characteristic SED in each pixel. The 

^ ! convolution kernels and the IDL packages used to construct and use them are 

^ ! made publicly available. 

\6 '■ 

O , Subject headings: Astrophysical Data, Data Analysis and Techniques, Astronom- 

I ! ical Techniques 

> ' 

^: 

■ 1. Introduction 

Spectral energy distribution studies of astronomical objects provide insight into the 
ongoing physical processes. In order to achieve a wide range of wavelengths, it is often 
necessary to combine observations from cameras with very different PSFs, in some cases 
with full width half maximum (FWHM) differing by factors of 100. Direct comparison 
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(e.g., ratios) of images with structures on multiple angular/spatial scales obtained with 
different PSFs can result in unphysical intensity ratios (i.e., colors). To preserve colors, 
intensity ratios should be calculated from images with a common PSF. We therefore require 
convolution kernels that will transform the images taken with several instruments into a 
common PSF, so we can generate image cubes (i.e., a collection of images expressed in the 
same sky coordinates grid) in which each pixel corresponds to the same sky region for all 
the cameras used. 

By "camera" we refer to the combination of telescope optics and physical detector, 
including the effect of atmospheric "seeing" if applicable. The PSF \l/j(x, y, x', y') of a camera 
i gives the measured intensity at (x, y) produced by a point source with unit flux at the point 
{x',y'), where we use the Cartesian coordinates {x,y) to denote positions in a small region 
of the sky. With this definition, the PSF has normalization 



for any source position {x',y'). 

It is often possible to approximate the PSF (denoted as \1/ from now on) as constant 
across the useful field of view of the camera, so '^j{x,y,x',y') = ^^(a; — x',y — y'). The 
observed image /^(x, y) will then be the convolution of the source S{x, y) with the PSF : 



Clearly, given two cameras A and 5, with (different) PSFs ^ a and '^Bi the images 
obtained of an astronomical object will be different, even if the spectral response of the 
cameras were identical. A convolution kernel is a tool that transforms the image observed 
by one camera into an image corresponding to the PSF of another camera. 

The convolution kernel K{A =^ B} from camera A to camera B should satisfy 



lB{x,y)= / lA{x\y')K{A^B}{x-x',y-y')dx'dy'={lA^K{A^B}){x,y), (3) 



where Ia and Is are the observed images by the cameras A and B respectively. 

The actual PSF of an instrument will show variations with the source color, variations 
along the field of view, changes over time, and deviations from rotational symmetry. Full 
2D characterization of a PSF is extremely challenging, and its extended wings are often not 
well determined. To take into account deviations of PSFs from rotational symmetry would 
require separate kernels K{A =^ B}{'-f) for each relative orientation 7 of cameras A and B. 




(1) 




x-x\y- y')dx'dy' = {S -k ^j) (x, y). 



(2) 




- 3 - 



In the present study, the actual PSFs are close enough of having rotational symmetry that 
such additional complication is not justified (see §3 for a detailed study). The current work 
assumes that "^j^ and can be approximated by rotationally symmetric functions. 



Using a different technique, lAlard fc Luptoru ( 1l998l ) presented a method for finding 



optimal kernels to convolve (to a common resolution) images of a sky region taken with 
a sin gle camera un d er diff erent seeing conditions. Using techniques similar to those used 



here, iGordon et al.l ( l2008l ) created kernels for the Infrared Array Camera (IRAQ and the 



Multiband Imaging Photometer for Spitzer (MIPS) cameras of Spitzer, and lSandstrom et al. 



fl2009l ) created kernels for the Infrared Spectrograph (IRS) of Spitzer. In the present work, 
we use the latest characterization of the PSFs of the cameras on the Spitzer Space Telescope 
(IRAC and MIPS), Herschel Space Observatory (Photocamera Array Camera and Spec- 
trometer for Herschel (PACS) and Spectral and Photometric Imaging Receiver for Herschel 
(SPIRE)), Galaxy Evolution Explorer (GALEX) (FUV, NUV), and the Wide-field Infrared 
Survey Explorer (WISE) (W1-W4). Additional PSFs (including those characterizing the IRS 
spectrograph on board Spitzer; the PACS spectrographs on board the Herschel Space Ob- 
servatory; and all-sky images produced by IRAS, AKARI, and Planck) are currently being 
constructed and will be added to the kernel library. We also include a family of analytical 
PSF profiles that are commonly used to model the PSFs for ground based telescopes. We 
construct the set of kernels to transform from different instrumental PSFs into Gaussian 
PSFs, another form that is widely used. We find an optimal Gaussian PSF for each camera: 
i.e., a Gaussian PSF that is sharp enough to capture the angular resolution of the camera, 
and wide enough to be robust against image artifacts. 

Additional image processing (i.e., co-adding image^ or different data reduction schemes 



i.e., the Scanamorphos pipeline ( lRousselll201ll ) for the Herschel images) will alter the PSF 



and thus new kernels should be constructed using the effective PSF. 

This article is organized as follows. In §2 we describe the generation of convolution 
kernels, in §3 we describe the PSF used, and in §4 we describe the kernel construction 
strategy. The performance of the kernels is examined in §5, and in §6 we describe a set of 
Gaussian kernels that are compatible with the different instruments. In §7 we describe the 
kernel usage and show their performance on NGC 1097, and we summarize in §8. 

All the kernels, IDL routines to use the kernels and IDL routines to make new kernels, 
along with detailed analysis ot tbe generated kernels, are publicly availableU 



^For example, for 2MASS co-added images the images undergo an additional smearing with a kernel 
whose size corresponds to a detector pixel. 

^See http://www.astro.princeton.edu/~draine/Kernels.html Kernels for additional cameras and updates 
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2. Convolution Kernels 

Given two cameras A and B, witli PSFs and we seek K{A =^ B} tliat fulfills 
equation ([3]) for any astronomical source. Thus, 

{S ^ ^b) = Ib = {I a * K{A B}) = {S*^A^ K{A B}) (4) 

for any astronomical source S, so the convolution kernel must satisfy 

<i/B = {^A^K{A^ B}). (5) 

With the normalization condition given by equation ([T]), the kernel must have 

K{A B}{x, y)dxdy = 1. (6) 

We can easily invert equation in Fourier space; taking the two-dimensional Fourier 
transform (FT) of equation ([5])|f| we obtain 

FT (^s) = FT {^A * K{A B}) = FT {^a) x FT {K{A B}) . (7) 
This can be inverted to obtain 

K{A ^B} = FT-' (fT(v1/^) x Jj^^ , (8) 

where FT and FT"' stand for the Fourier transform and its inverse transformation respec- 
tively. 

Equation (jH]) provides a condition for the existence of such a kernel and a practical way 
of computing it. We can see that a condition for the existence of a kernel is that the Fourier 
components for which FT(\E'^) = should satisfy FT(\E's) = 0. Informally speaking, this 
means that the PSF of camera A must be narrower than the PSF of camera B. 

For each camera A, we identify a high-frequency cutoff kH,A as the highest spatial 
frequency for which FT(\E^yi) is appreciable by setting 

FT{^A){kH,A) = 5 X 10-=^ X max[FT(^^)]. (9) 



will be included when new PSF characterizations become available. 

■^Computing the two-dimensional Fourier transform of rotationally symmetric functions is mathemati- 
cally equivalent to making a decomposition in the one-parameter family of Bessel functions. However, the 
existence of the fast fourier transform (FFT) algorithm makes it numerically more efficient to perform the 
decomposition in the family of Fourier modes. 
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The cutoff frequency kn^A can be normalized as kn.A = i^a^ 27r/FWHMA. The values of ka 
are in the range of 1.08 - 1.46 and are given in Tabled! 

Equation ([HD also provides insight of a possible problem in computing kernels. PSF 
Fourier transforms do not have significant power at (spatial) frequencies above the 27r/FWHM. 
The high-frequency components of the FT will be small, introducing large uncertainties when 
inverted. 

A possible way of avoiding problems in the high frequency components of the kernel is 
to introduce a filter Ja in the kernel construction: 



K{A ^ B} = FT~^ FT{mB 



1 



(10) 



where fA is a suitable low-pass filter. Because this differs from equation ([H]), it is clear that 
a kernel satisfying equation ( fTOl) will, in general, not be an exact solution to equation (E]). 
However, we can expect that if the filter does not remove significant power from the high 
frequency components of either FT{^>a) or FT(\I/b), then the kernel computed will be a 
good approximate solution to equation ([8]). 



We use a filter Ja of the form 



fA{k) 



1 

1 

2 





X 



1 + cos ( TT X 



k-k 



L,A 



kH,A-kL,A 



for k < kL,A 
for kL,A <k< kH,A, 
for kn A k 



where we set k^^A = 0.7x /c^^^^. Note that the cutoff frequency of our filter /a depends only on 
camera A, since small values of FT(\I/b) have no negative impact in the kernel construction. 
In principle, any smoothly varying function that is close to 1 in the low-frequency range and 
goes to zero in frequencies larger than kH,A should work as well. We have experimented using 
several smoothing functions and find that the particular form of fA{k) given by equation 
( ITTj) . with these particular choices of k^^A and ku^Ai gives excellent results. More details of 
the smoothing function effects can be found in §5. 

If the cutoff frequencies associated with cameras A and B are such that kfj B ^ ^h^a x 0.7, 
then the filter has little effect and the resulting kernel from equation ffTOj) will satisfy 
equation (|3]) to a very good approximation. In the regime kn^A x 0.7 < kn^B < kH,A the 
kernel K{A =^ B} will depend on the exact form of the filter Ja used. Most of the kernels 
in this regime have good performance; the performance of such "filtered" kernels will be 
evaluated in §5. A limiting case is a kernel K{A =^ A} that transforms a PSF into itself; it 
is the Fourier transform of the filter Ja- 
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In the cases of kH,A < kH,B (convolving into narrower PSFs), use of the filter Ja allows 
one to compute convolution kernels, but their performance can be poor. This will be further 
discussed subsequently. 



3. Instrumental point-spread functions 



We generate appropriate kernels for the measured PSFs of the cameras on Spitzer, 
Herschel, GALEX, and WISE and for certain analytical PSFs that are used in the literature: 



IRAC- The Infrared Array Camera (IFazio et al.ll2004l ) has four infrared bands, centered 
at 3.6/im , S.O/xm , 5.8/im , and 8.0/im . We use the in-flight extended point response functions 
(PRF)0 The core portion of the PSF was made with 300 observations of a calibration source, 
with different exposure times, combined into a high dynamic range image. Observations of 
the stars Vega, eEridani, Fomalhaut, eindi and Sirius were used in the construction of the 
extended wings of the PSfI. 



MIPS.- The Multiband Imaging Photometer for Spitzer (IRieke et al.l 120041) has three 
photometric infrared bands, c entered at 24/im , 7 /xm , and 160/xm . Following lEngelbracht et al 
( 120071 ): iGordon et all ( 120071 ): IStansberry et al.l (120071 ). we generate theoretical PSFs for the 
MIPS cameras, assuming a blackbody source at T = 25K. The PSFs are generated using the 
software STinyTim0 in a 0.5" grid, and they are smoothed with a square kernel with sides of 
4.5" , 13.5" , and 25.5" for the bands at 70yum , 100/xm , and 160/xm respectively. The smooth- 
ing sizes correspond to 1.6, 1.35, and 1.8 times the camera detector pixel size, and they 
should cause the core of the theoretical PSF to be in close agreement with the calibration 
point- source images. 



PACS.- The Photocamera Array Camera and Spectrometer for Herschel ( jPoglitsch et al. 



2010l ) has three photo metric infrared bands, centered at 70/xm, lOO/imand IGOum . We 
use the in-flight PSF dCeis fc Lutjboioi Ihutzl boioi iMuller fc The PACS icdboiol ): the 
core was defined by observations of the star a Tau and the asteroid Vesta, with extended 
wings reconstructed from (saturated) observations of Mars, Neptune, IK Tau, and the Red 
Rectangle. The azimuthally averaged encircled energy fraction of the PSFs out to 1000" were 
obtained from the HCSS/HIPE software. 



"'in the current work we always use the full response of the optical systems including the camera effects, 
and for simplicity we will not make any further distinction between PSFs and PRFs 

^The PSFs available at: |http://ssc. spitzer. caltech.edu/ 



^STinyTim is available at: 



http://irsa.ipac.caltech.edu/data/SPITZER/docs/dataanalysistools/tools/contributed/genera] 
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^SPIRE.- The Spectral and Photometric Imaging Receiver for Herschel ( iGriffin et al. 

2010 ) has three photometric far-infrared bands, centered at 250/im, 350/xmand 5OO/im0 



GALEX.- The Galaxy Evolution Explorer (IMartin et al. has two ultraviolet bands, 
FUV (1350 - 1750 A) and NUV (1750 - 2800 A)@ 



WISE.- The Wide-field Infrared Survey Explorer (jWright et al.ll2010l ) has four photo- 
metric infrared bands, centered at 3.4/im , 4.6/xm , 12/im , and 22/im . The PSF shape varies 
significantly over the focal plane due to distortion from the telescope optics; so a library of 
PSFs corresponding to a 9 x 9 grid of locations on the focal plane was determined. For each 
camera, we average the 81 different PSFs to generate a single PSF per camera. 

As a way of measuring the departure of a PSF from rotational symmetry, we define 
an asymmetry parameter: 



gj = j j \ - C[^j\ I dxdy, 



(12) 



where C[\E'j] is the azimuthally- averaged version of ^ j (i.e., is a PSF with rotational 

symmetry and the same radial profile as j)- The PSFs used have g < 0.2. 

In Table [1] we have a summary of the properties of the different PSFs: the camera 
Rayleigh diffraction angle, the PSF FWHM, the radius containing 99% of the PSF energy, 
the (normalized) high-frequency cutoff ka used in the filter and the anisotropy parameter 
g. The PSF radial profiles [out to '^{9) ~ 10~^\l/(0)], and enclosed power are plotted in 
Figures [1] and [21 

For each PSF the radii containing 25%, 50%, 65%, 80%, 90%, 95%, 98%, 99%, 99.5%, 
and 99.9% of the total power are given in Table O Table [3] gives the enclosed power for 
selected radii. 



^We use the in-flight 1.0" PSF maps from "ftp: / / ftp.sciops.esa.int /pub/hsc-cahbration / SPIRE / 



In-flight PSFs arc available at http://www.galex.caltech.edu/researcher/techdoc-ch5.html 
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Table 1: Basic Instrument Information 





Raylcigh diff. 


Measured 


99% of energy 




Asymmetry 


Camera 


limit" (" ) 


FWHM (" ) 


radius (" ) 








1 n4 




fi9 ^0 


1 9Q 
X . z y 


n 1 fi 

U. lU 


TR AC 4 5//m 


1.31 


1.81 


64.46 


1.26 


0.17 


TR AC 5 8//m 


1.68 


2.11 


133.55 


1.20 


0.19 


IRAC S.O^m 


2.30 


2.82 


114.20 


1.19 


0.18 


MIPS 94inn 

iVXXX O ^^LLLLL 


fi Q3 


6 43 


224.53 


1 05 


08 


MIPS 7()iim 

XVXXX 1 UlAlJ-l 


20.90 


18.74 


461.44 


1.12 


0.05 


MIPS IfiO/zm 

XVXXX k_J XVJVJtXlll 


45.62 


38.78 


678.77 


1.10 


0.05 


PACS 7n/;iTi 


5.11 


5.67 


249.81 


1.23 


0.20 


PACS 1 nn//m 




7 04 


350 63 


1 19 

X . X C/ 


90 


PACS 1 fiflz/m 


1 1 70 

X X . ( u 


11.18 


41 7 36 


1.21 


90 


OJr IXVJZj ZtJUi-Alll 




1 8 1 fi 

lo. lo 


90^ 07 


1 1 f! 

1 . lU 


u. ly 


O-T XXL-Lj 00\J fJjilL 


9"^ 


94 88 


1 Q9 47 


1 1 ^ 


U. xO 


LJJT XX\jX1j 0\J\J Lllll 


99 


"ifi riQ 

OU.U J 




1 1 fi 


vj. X y 


CAT,F,X FT TV 


n 08 


4.48 


^0 98 


1.26 


07 


r" AT FY NTTV 


nil 

U. 11 


J.UO 


oy .ou 


1 '^9 
1 .oz 


0^ 


VV lO-Cj O.OO^IIl 


9 11 
Z. 1 1 


^ 7Q 


1 Q 1 n 


1 90 
l.ZU 


17 


AA^TSF A Rn//m 
VV lO-Cj ^l.OU^IIl 


9 SQ 


O.O ( 


ly.uo 


1 11 


1 


VV iO-Cj ii.OD/im 


7 97 


D.DU 


ly.OD 


1 91 
i.Zo 


19 
U. iZ 


WTSF 22 1 /;m 

V V XkJXJ Zj . X Lii 111 


13.90 


11.89 


35.15 


1.05 


0.07 


PpncG 1 9" 




1 9 nn 


1 ^ 41 


1 11 




u.u 


Pcmcc on" 
V^cIIIdd ZU 




9n nn 
zu.uu 


ZO.Do 


1 11 




u.u 






9'? nn 
Zo. uu 


90 "^"^ 


1 11 
i.oo 




u.u 






98 nn 
Zo. uu 


oo.yo 


1 11 

1.00 




u.u 


Gauss 40" 




40.00 


51.33 


1.33 


0.0 


Gauss 50" 




50.00 


64.05 


1.33 


0.0 


BiGauss 0.5" 




0.50 


0.90 


1.37 


0.0 


BiGauss 1.0" 




1.00 


1.79 


1.37 


0.0 


BiGauss 1.5" 




1.50 


2.69 


1.37 


0.0 


BiGauss 2.0" 




2.00 


3.57 


1.37 


0.0 


BiGauss 2.5" 




2.50 


4.44 


1.37 


0.0 


Moffat 0.5" 




0.50 


1.39 


1.46 


0.0 


Moffat 1.0" 




1.00 


2.77 


1.46 


0.0 


Moffat 1.5" 




1.50 


4.12 


1.46 


0.0 


Moffat 2.0" 




2.00 


5.45 


1.46 


0.0 


Moffat 2.5" 




2.50 


6.74 


1.46 


0.0 



" We take the Rayleigh diffraction limit as 1.22 x [central A] / [telescope diameter]. 
We define k as {kn x FWHMyi)/(27r) where kn is the high-frequency cutoff (see text for details). 
The parameter g is a measure of the departure of a PSF from rotational symmetry [see ea. (|12l) ]. 
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Table 2: PSF Enclosed Energy Percent (%) in Selected Circular Apertures 
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100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


BiCauss 1.0" 


99.9 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


BiGauss 1.5" 


98.6 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


BiGauss 2.0" 


95.7 


99.9 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


BiGauss 2.5" 


89.9 


99.5 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


Moffat 0.5" 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


Moffat 1.0" 


98.7 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


Moffat 1.5" 


96.1 


99.5 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


Moffat 2.0" 


90.9 


98.7 


99.7 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


Moffat 2.5" 


83.0 


97.7 


99.3 


99.8 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 


100 
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Table 3: Radii (" ) Enclosing Selected Percents of the Total Power 



Camera 



IRAC 3.6^im 
IRAC i.b^im 
IRAC 5.8fim 
IRAC 8.0/^m 



MIPS 24^m 
MIPS lOi^im 
MIPS leO^im 



PACS 70^m 
PACS 100/im 
PACS 160/im 



SPIRE 25Qfim 
SPIRE 350nm 
SPIRE SOO^im 



GALEX FUV 
GALEX NUV 



WISE 3.35^m 
WISE i.eO^im 
WISE ILSe^im 
WISE 22.1^m 



Gauss 12" 
Gauss 20" 
Gauss 23" 
Gauss 28" 
Gauss 40" 
Gauss 50" 



BiGauss 0.5" 
BiGauss 1.0" 
BiGauss 1.5" 
BiGauss 2.0" 
BiGauss 2.5" 



Moffat 0.5" 
Moffat 1.0" 
Moffat 1.5" 
Moffat 2.0" 
Moffat 2.5" 



25 



50 



65 



80 



Percent 

(%) 
90 95 



98 99 99.5 99.9 



0.74 
0.75 
0.99 
1.16 



1.19 
1.24 
1.97 
2.21 



1.60 
1.73 
3.17 
3.32 



2.44 
2.68 
16.0 
6.70 



5.83 
7.06 
42.3 
18.9 



12.9 
14.7 
75.3 
39.5 



36.0 
37.9 
114 
82.9 



62.5 
64.5 
134 
114 



88.1 
90.8 
145 
134 



2.73 
8.02 
16.3 



4.93 
14.4 
29.0 



9.18 
25.5 
52.8 



12.2 
35.6 
71.9 



29.4 
80.2 
155 



60.4 

158 
287 



130 
318 
518 



225 
461 
679 



368 
628 
802 



2.41 
2.83 
4.50 



4.23 
4.87 
7.66 



6.47 
7.50 
11.3 



12.5 
14.1 
20.8 



37.5 
43.5 
53.3 



85.3 
105 
133 



165 
227 
294 



250 
351 
417 



378 
477 
524 



6.39 
8.49 
12.4 



10.2 
13.4 
19.5 



13.1 
17.0 
24.6 



20.8 
23.7 
33.5 



36.6 
43.9 
60.3 



66.1 
75.0 
85.9 



138 
122 
137 



205 
192 
198 



382 
397 
411 



1.59 
1.82 



2.49 
2.90 



3.17 
3.71 



4.27 
5.08 



7.05 
9.34 



19.0 
20.9 



38.3 
33.9 



50.3 
39.6 



59.7 
43.4 



2.03 
2.29 
3.02 
5.01 



3.26 
3.69 
5.48 
8.87 



4.17 
4.75 
7.66 
14.6 



5.67 
6.36 
10.2 
19.1 



8.36 
9.00 
12.8 
22.3 



12.5 
13.0 
15.1 
25.4 



17.2 
17.2 
18.1 
31.2 



19.1 
19.1 
19.6 
35.1 



20.2 
20.2 
20.5 
38.5 



3.86 
6.44 
7.41 
9.02 
12.9 
16.1 



6.00 
10.00 
11.5 
14.0 
20.0 
25.0 



7.38 
12.3 
14.1 
17.2 
24.6 
30.8 



9.14 
15.2 
17.5 
21.3 
30.5 
38.1 



10.9 
18.2 
20.9 
25.5 
36.4 
45.5 



12.5 
20.8 
23.9 
29.1 
41.5 
51.9 



14.2 
23.7 
27.3 
33.2 
47.4 
59.2 



15.4 
25.7 
29.5 
36.0 
51.3 
64.1 



16.5 
27.5 
31.6 
38.5 
54.9 
68.4 



0.17 
0.33 
0.50 
0.66 
0.83 



0.26 
0.52 
0.78 
1.04 
1.30 



0.32 
0.64 
0.97 
1.29 
1.61 



0.41 
0.81 
1.22 
1.62 
2.03 



0.50 
1.00 
1.50 
2.00 
2.50 



0.60 
1.20 
1.80 
2.40 
3.00 



0.76 
1.53 
2.29 
3.04 
3.79 



0.90 
1.79 
2.69 
3.57 
4.44 



1.02 
2.04 
3.05 
4.04 
5.01 



0.18 
0.36 
0.54 
0.71 
0.89 



0.29 
0.57 
0.86 
1.14 
1.43 



0.36 
0.72 
1.08 
1.44 
1.81 



0.47 
0.94 
1.41 
1.88 
2.34 



0.60 
1.21 
1.81 
2.42 
3.02 



0.76 
1.53 
2.29 
3.05 
3.80 



1.07 
2.14 
3.19 
4.23 
5.27 



1.39 
2.77 
4.12 
5.45 
6.74 



1.72 
3.43 
5.09 
6.69 
8.29 



132 
134 
154 
152 



738 
894 
950 



712 
711 
712 



499 
511 



75.1 
58.7 



21.4 
21.4 
21.5 
42.1 



18.5 
30.9 
35.5 
43.2 
61.5 
76.3 



1.25 
2.50 
3.71 
4.86 
5.96 



2.26 
4.50 
6.70 
8.90 
11.1 
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Fig. 1. — PSF radial profiles for the 20 cameras considered here (see text). All PSFs are 
shown out to ^(^) ^ 10"^\l/(0), with the exception of WISE (W1-W4), for which we have 
data only down to ^ 0.002^(0) 
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Fig. 2. — Fraction of the power enclosed by a circle of radius Q for the 20 cameras considered 
(see text). 
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For convenience, we add several families of PSFs that are often used for ground-based 
optical and radio telescopes. For each analytical profile, we generate the kernels for a range 
of FWHM values. We consider the following analytical profiles: 



Gaussians.- We use Gaussian PSF of the form 

1 



2 vr 0"' 



exp 



2^2 



(13) 



where the FWHM = 2x a a/2 In 2. We generate kernels with 5" < FWHM<65'' 



Optical (SDSS).- In the SDSS survey, it is found that a good approximation to the 
telescope PSF is given by the sum of two Gaussians. The two components have relative 
weights of 0.9 and 0.1, and the FWHM of the second component is twice that of the first. 
We use a family of PSFs of the form: 

1 / a2 



^{6) = 0.9 



2 TT 0"' 



_02 

expi ^ 1+0.1 



2 7r(2(j)2 



exp 



2(2a) 



(14) 



where the FWHM = 1.01354 x 2 x a ^2 In 2. We generate kernels with FWHM = 0.5, 1.0, 
1.5, 2.0, and 2.5". 



Optical (General) .- iMoffatI ( 1l969l ) proposed PSFs of the form: 



(/3_l)(2i/^_l) 



1 + (2 



1 -P 



(15) 



where /3 is a parameter, and FWHM = 2 6'o, Following iRacind (119961 ) we use PSFs of the 
form: 

m[r) = 0.8 X Mjie) + 0.2 X M2(^), (16) 

where the same is used in My and M2. We generate kernels with FWHM = 29q = 0.5, 
1.0, 1.5, 2.0, and 2.5". 



4. Kernel Generation 



Generation of the convolution kernels was accomplished as follows: 



4.1. Input the PSF and correct for missing data. 



When an input PSF image has missing data pixels (like the current SPIRE PSFs), we 
iteratively estimate the missing values. 
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We start by replacing the missing data pixels by a value of in the original image. 
We compute a smoothed image by convolving it with a (normalized) Gaussian kernel oc 
exp[— (6'/26'o)^], with Oq equal to 2 pixels. We replace the original image missing data pixels 
by the value they have in the convolved image (the original data are not altered). We iterate 
the convolution and replacement steps 5 times. The resulting image has all the missing data 
points replaced by a smooth interpolation. This technique produces robust results, even if 
we have missing data in a patch of a few contiguous pixels. 



4.2. Resample the PSFs. 

Each PSF comes in a grid of different pixel size. We transform each PSF into a grid of 
a common pixel size of 0.20" using the IDL procedure congrid, using the cubic convolution 
interpolation method with a parameter of -0.5. The 0.20" pixel size capture all the details 
on the instrumental and Gaussian PSFs. We also pad with the resulting images into an 
odd-sized square array if needed. 

We use a grid of a common pixel size of 0.10" to regenerate the kernels from optical 
PSFs into IRAC cameras. 



4.3. Center the PSFs. 

To determine the image center, we smooth the image with a 5 pixel radius circular 
kernel, and locate the image maximum. In some PSFs the maximum value is achieved over 
a (small) ring. To avoid possible misidentification of the real image center, we take the 
centroid of all the pixels that satisfy: 

max[\l/] - \l/(x,?/) 4 



max[\E'] 



< 5 X 10-^ (17) 



4.4. Circularize the PSFs. 

In order to make a rotationally symmetric PSF, we average over 2^^ rotations of the 
image through angles /3„ = x 360°/2^^ for n = 1,2,3, ...,2^^ = 16,384, producing a PSF 
image that is invariant under rotations of any angle that is a multiple of 360°/2^^ = 0.022° 
(i.e., is as rotationally symmetric as we can numerically expect). 

Computing 2^^ rotations naively would be extremely time-consuming, but the final result 
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can be in fact computed performing only 14 rotations, as follows. 

We start by rotating \[^by an angle 6i = 360°/2^ = 180°, producing an image -Ri, and 
computing their average Ri = ^ x + Ri] . Clearly, Ri is now invariant under rotations of 
= 36072^ = 180°. 

We continue this procedure iteratively. We rotate -Ri by an angle 62 = 360°/2^ = 90°, 
producing an image R2, and set R2 = ^ x [Ri + R2]. R2 is invariant under rotations of 
9i = 360°/2^ = 180° and 62 = 360°/2^ = 90°; i.e. it is invariant under rotations of any angle 
that is a multiple of 360°/22 = 90°. 

We iterate this procedure 14 times; the last average calculated image R14 is the desired 
rotationally symmetric PSF. 

We further set to all the pixels that lie outside the largest circle included in the square 
image, since those regions would correspond to areas with partial image coverage. If there 
are pixels with (very small) negative values (due to the noise in the original PSFs) we set 
them to 0. 

In order to have a more stable algorithm, the previous rotations are performed in reverse 
order (i.e., starting with the smallest angles). 

All the remaining steps in the kernel construction should preserve the rotational sym- 
metry in the images. A way of estimating the noise induced by some steps (e.g., computing 
Fourier transform) is to compute the departure from rotational symmetry in the resulting 
image. Circularizing helps to reduce numerical noise, and will be performed after every 
step that can potentially decrease the image quality. When a circularization is performed 
to a rotationally symmetric image, the asymmetry parameter g of the resulting image is 
smaller than 0.0008 (i.e., the circularization procedure itself induces very small departures 
from rotational symmetry). 

4.5. Resize the PSF. 

We trim (or pad with 0) all the PSFs into a common grid, to be able to compute all the 
convolution kernels using only one Fourier transform per PSF. We choose a grid size that is 
large enough to contain most of the power in each PSF. We also optimize its size to make 
the FFT algorithm as efficient as possible (minimal sum of prime factors) . The adopted grid 
size is 729" x 729" , giving an image size of 3645 x 3645 pixels. 
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4.6. Compute the Fourier transform of the PSF FT(\1/). 

We use an efficient FFT algorithm. Since the PSFs are invariant under reflections, 
X -H- —X, their Fourier transform should be real. We impose this real condition to reduce 
the numerical noise introduced by the double-precision FFT algorithm. 



4.7. Circularize the FT(^). 

Using the procedure as before, we circularize the FT(\E'). In principle, they should 
already be rotationally symmetric, but numerical noise in the FFT algorithm makes them 
slightly non rotationally symmetric. 



4.8. Filter the FT(^). 



We filter the highest frequencies in each FT(\i/). We use a filter of the form 



m 



exp 




1.8249 X 



for k <ka 
for ka < k < ki3, 
for kp < k 



where we set ka = 0.9 x kjj. The factor 1.8249 is chosen so that 0(0.5 x [kp + ka)) = 0.5. For 
each camera we choose = 4 x (27r/FWHM). We tested several filter functions and found 
that the particular form given by equation ( 1181) provided the best results. Each PSF has 
structure at spatial wavelengths comparable with the FWHM, so the Fourier components 
with frequencies much higher than this cannot be important. The Fourier components 
removed by this filter were mainly introduced by the original image resizing algorithm. In 
the following discussion, we let FT^(\E') = x FT(\E'). 



4.9. Invert the FT(^). 



We evaluate 1/FT^{^) at the points where FT^{<iJ) ^ 0. We set 1/FT<^(^) = at the 
remaining points (which will be filtered soon). 
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4.10. Compute the FT of the filtered kernel. 

We compute the FT of the filtered kernel using the filter Ja from equation ffTTl) : 

FT{K{A ^ B}) = FT^i^B) x JjT^ (19) 
for all the appropriate combinations {A,B). 



4.11. Compute the kernels. 

We compute the inverse Fourier transform to the previously calculated FT{K). We 
again impose the condition that K must be real. 



4.12. Circularize the kernels. 

Using the procedure as before, we circularize the kernels. Again, they should already be 
rotationally symmetric. Numerical noise in the inverse FFT algorithm makes them slightly 
asymmetric, but this is easily corrected. 



4.13. Resample the kernels. 

All the computed kernels are given in a grid of a common pixel size of 0.20" , but will 
be needed in grids of different pixel sizes. Again, we resample the kernels using the IDL 
procedure congrid, using the cubic convolution interpolation method with a parameter of 
-0.5. 



4.14. Final trim of the kernels. 

We trim each kernel to a smaller size (to speed up further convolution) such that it 
contains 99.9% of the total kernel energy. Moreover, we use a square grid with an odd 
number of pixels so that the kernel peaks in a single central pixel. 
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4.15. Circularize the final Kernels. 

We finish the kernels by circularizing them again, to remove the small noise introduced 
in the resampling process. 

All previous calculations were done in double precision to reduce numerical noise. 



5. Kernel Performance 

For each generated kernel, we compute * K{A =^ B} the convolution of \E'a and 
K{A =^ B}, and compare it with ^bE For a perfect kernel, both quantities should coincide 
at all radii. 

Figure [3] shows the analysis of K{M2A =^ S25O}0 This kernel shows good behavior: it 
transforms from a camera with FWHMm24 = 6.5" into a camera with FWHMs25o = 18.2". 
This kernel essentially spreads the energy of the core of MIPS 24/imPSF into a larger area. 
The filter /m24 has no effect on the construction of this kernel, because FT(\E'b) ~ for 
k > 0.7 X kff^A 

The left panel of Figure |3] compares the integrated power of the PSFs. It includes '^M2i 
(dot-dashed line), \E's25o (solid line), and \E'M24*-f^{M24 =^ S250} (dashed line). For an ideal 
kernel the solid line and the dashed line should coincide. The departures in this case are 
very small. 

The power per unit radius is proportional to ^\l/(6'), or OK {9). The right panel of Figure 
|3] shows 9'^ and 9K (normalized to the maximum value). The lower part of the right panel 
includes four traces: 9 x \1/m24 (dot-dash line), 9 x \E's25o (sohd line), 9 x (\E'm24 * -f^{M24 =^ 
S250}) (dashed line) and 9 x _ft'{M24 ^ S250} (dotted line) for visualization of the kernel 
behavior. For an ideal kernel the solid line and the dashed line should coincide. The upper 
part of the right panel is a plot of the difference between 9 x ^^3250 and 9 x (\E'M24*-f^{M24 ^ 
S250}). For an ideal kernel this graph should be 0. For this example (i^{M24 ^ S250}), 
the ^M24* A'{M24 ^ S250} reproduces the SPIRE 250/imPSF to within 0.01 x ^s25o(0). 

In Figures Hand [5] we analyze the kernels Jr{M70 S250} and i^{S250 M70}. Their 



^It can be easily shown that the radial profile of {'^ a*K{A ^ B}) is the same whether ^ a is rotationally 
symmetric or not, so for simplicity, we use the circularized version of a in the comparison. 

10 In our PSF and kernel notation, we will abbreviate I, M, P, S, GAL, and W for IRAC, MIPS, PACS, 
SPIRE, GALEX, and WISE respectively, and will omit the /nm symbol. 
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Fig. 3. — Performance of the filtered kernel K{M24 =^ S250}. W- is the integral of the 
negative values of the kernel, and D is the integral of the absolute value of the difference 
between the target PSF and the PSF reproduced by the kernel [see eq. f l20|) ]. See the 
electronic edition of the PASP for a color version of this figure. 



construction is more challenging since both cameras have similar FWHM: FWHMmto = 
18.7" and FWHMg25o = 18.2". These kernels have to redistribute the energy within the 
core and Airy rings of the PSFs. The plotted quantities are similar to those in the right 
panel of Figure [3] for K{M24 =^ S250}. In this case, the kernels still perform well, but they 
have large areas with negative values. 

One measure of kernel performance is its accuracy in redistribution of PSF power. We 
define 

D = J J I - ^ fi} ^ I rfxdy. (20) 

A kernel with perfect performance will have D = 0, and normalization of the PSFs requires 
D < 2. D measures how much flux has not been redistributed correctly. Good kernels 
have small D values: D{K{M24: S250}) = 0.011. In Table H we give D for the kernels 
constructed. 

A second quantitative measure of kernel performance is obtained by studying its negative 
values. We define 

W± = ^ f I {\ K{A =^ B} I ±K{A B}) dxdy. (21) 
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Fig. 4. — Performance of the kernel K{M70 =^ S250}. Because of the large value W- = 2.14, 
convolution MIPS70/im^ SPIRE250/imis not recommended. See the electronic edition of 
the PASP for a color version of this figure. 



Flux conservation requires that Wj^ = 1 + W^. In general, kernels will have > 0. 
Well-behaved kernels have small W- values: W-{K{M24: S250}) = 0.07. The integral of 
I K{A =^ B} I is [1 + 212r_], so a kernel with a large value of W- could potentially amplify 
image artifacts. Additionally, a kernel with large can generate areas of negative flux 
near point sources if nonlinearities are present, or if the background levels were subtracted 
incorrectly. Table [5] lists the IV_ values for the kernels constructed. The values in Table 
|5] were computed numerically. Due to finite grid resolution the numerical values may be 
off by a few percent in some cases. This can be seen from the IV_ values computed for the 
self- kernels. The self-kernels are simply the Fourier transform of the filter and should 
therefore be the same (1.15) in all cases. However, in Table [5] the W- for the smallest PSFs 
are larger than 1.15 by as much as 0.03 (e.g., 1.18 for PACS 70/im). 

Essentially, there are two sources of oscillations due to the filter and the need 
to remove energy from some region to relocate to another region (when the target ^E'b is 
narrower than or has less energy in some annuli). 

The kernels between cameras with similar FWHM also have oscillations. Using a softer 
filter J'a would reduce the oscillatory behavior, giving smaller values of W^, but would also 
produce worse matched PSFs (larger values of D). The particular form of filter J'a used in 
the present work is a good compromise between having good PSF matching and moderate 
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20 40 60 80 100 120 140 
9 (arcsec) 



Fig. 5. — Performance of the kernel ii'{S250 =^ M70}. With W- = 0.61, this kernel is safe 
to use. \l/s25o*-ft'{S250 =^ M70} deviates from \1/m7o less than 2%. See the electronic edition 
of the PASP for a color version of this figure. 



oscillatory behavior. 

In Figure [5] we analyze a kernel of particular interest: ii'lMlGO =^ S500}. Despite 
both cameras having similar FWHM (FWHMmi6o/FWHMs5oo = 1-07), their extended wings 
are very different. The kernel A'{M160 =^ S500} is particularly badly behaved, with 
large negative excursions, having = 2.81 and D = 0.17. In a convolution of MIPS 
160/im images of NGC 6946, some bright point sources generated regions with negative flux 
around them. We do not recommend using the kernel i^{M160 =^ S500}: if MIPS160/imand 
SPIRE500/im images need to be combined, we recommend using the PSF of MIPS 160 fim ( 
ii-jSSOO M160} has ^ 0.47 and D ^ 0.042) or some Gaussian PSF compatible with 
MIPS 160 /im, such as a Gaussian with FWHM=64" (see §6). 

Finally, in Figure [7] we analyze the kernel i^{M70 =^ M70}. This kernel is essentially 
FT(/m7o), and it illustrates the effect of all the kernel construction steps. All the kernels of 
the form K{A =^ A} are scaled versions of i^'{M70 =^ M70}, aside from small differences 
due to finite grids. All of the K{A =^ A} kernels have ~ 1.15 and D ^ 0.06. 

Even though kernels K{A =^ A} have = 1.15 and = 2.15, they do not amplify 
the noise that arises from astronomical sources. The image of an astronomical point source 
will be When we convolve the camera with a kernel of the form K{A =^ A}, the image 
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Table 4: D Value for Constructed Kernels 



From 


To 


M M M 
24 70 160 


P P P 
70 100 160 


S S S 
250 350 500 


GAL GAL 
FUV NUV 


w w w w 

3.4 4.6 11.6 22.1 


io.o 
14.5 
15.8 
io.U 


0.065 0.029 0.037 
0.065 0.029 0.037 
0.064 0.031 0.048 
0.064 0.030 0.047 


0.051 0.047 0.031 
0.051 0.047 0.031 
0.047 0.045 0.031 
0.047 0.046 0.032 


0.025 0.022 0.018 
0.025 0.022 0.018 
0.029 0.027 0.018 
0.026 0.024 0.018 


0.006 0.005 
0.006 0.005 
0.012 0.010 
0.012 0.007 


0.005 0.005 0.007 0.008 
0.005 0.005 0.007 0.008 
0.012 0.012 0.013 0.012 
0.007 0.007 0.009 0.009 


M24 

M ( U 

M160 


0.091 0.018 0.048 
NC 0.055 0.038 
NC NC 0.064 


0.213 0.068 0.015 
NC NC NC 
NC NC NC 


0.011 0.013 0.017 
0.102 0.015 0.014 
NC NC 0.169 


NC 0.352 
NC NC 
NC NC 


0.225 0.122 0.098 0.027 
NC NC NC NC 
NC NC NC NC 


P70 
PlOO 
P160 


0.006 0.019 0.049 
U.UiZ U.U14 U.U40 
NC 0.009 0.040 


0.045 0.010 0.014 
U.ioo U.UoU U.UiU 
NC NC 0.044 


0.011 0.014 0.017 
U.Uii U.U14 U.Uio 
0.013 0.014 0.015 


0.227 0.113 

IN i\y_^ 

NC NC 


0.064 0.045 0.027 0.027 
U.zUz U.iii U.Uo / U.Uzo 
NC NC NC 0.032 


S250 
S350 
S500 


NC 0.037 0.045 
NC 0.300 0.043 
NC NC 0.042 


NC NC NC 
NC NC NC 
NC NC NC 


0.058 0.009 0.014 
NC 0.060 0.012 
NC NC 0.055 


NC NC 
NC NC 
NC NC 


NC NC NC NC 
NC NC NC NC 
NC NC NC NC 


FUV 

NUV 


0.065 0.030 0.036 
0.065 0.030 0.037 


0.049 0.046 0.032 
0.051 0.046 0.032 


0.025 0.024 0.018 
0.027 0.024 0.018 


0.046 0.017 
0.071 0.028 


0.007 0.006 0.007 0.007 
0.008 0.006 0.006 0.007 


W3.4 
W4.6 
W11.6 
W22.1 


0.066 0.029 0.037 
0.066 0.029 0.037 
0.074 0.028 0.037 
NC 0.028 0.036 


0.080 0.048 0.032 
0.079 0.047 0.032 
0.117 0.056 0.031 
NC NC 0.163 


0.026 0.023 0.018 
0.026 0.023 0.018 
0.026 0.023 0.018 
0.028 0.023 0.018 


0.252 0.122 
NC 0.118 
NC 0.206 

NC NC 


0.053 0.034 0.017 0.007 
0.051 0.033 0.016 0.006 
0.109 0.063 0.040 0.005 
NC NC NC 0.094 



Notes." D is the integral of the absolute value of the difference between 

a target PSF and the PSF reproduced by the kernel [see eq. ([20| ]. 

We are abbreviating I, M, P, S, GAL, and W 

for IRAC, MIPS, PACS, SPIRE, GALEX, and WISE respectively. 

NC stands for not computed: the kernel performance would be too poor. 



of a point source will still be very close to since D{K{A =^ A}) ~ 0.06. This implies 
that the noise coming from a field of unresolved astronomical background sources will not 
change significantly when we convolve the image with K{A =^ A}. To verify this reasoning, 
we generate an image 5* having independent Gaussian noise in each pixel in a very fine 
(0.2" ) grid. We convolve S with \1/m7o to have a simulated observed image of the noise: 
O = ^M70- We further convolve O with K{M70 M70}: C = 0^ K{M70 =^ M70}. 
We found that | 1 — (C^)/(0^) |~ 10~^. While uncorrelated noise is not amplified by a 
self- kernel K{A ^ A}, imprecise characterization of the PSFs and camera artifacts can be 
amplified by kernels with large W- values. 

There is no single number that captures all of the characteristics of a convolution kernel, 
but we find that W- serves as a good figure of merit. Based on experimentation with various 
kernels, we regard kernels with < 0.3 to be very safe. Kernels with ~ 0.5 also appear 
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Table 5: Value for Constructed Kernels 



From 


To 


M M M 
24 70 160 


P P P 
70 100 160 


S S S 
250 350 500 


GAL GAL 
FUV NUV 


w w w w 

3.4 4.6 11.6 22.1 


13.6 
14.5 
15.8 
18.0 


0.00 0.00 0.00 
0.01 0.00 0.00 
0.15 0.03 0.02 
0.11 0.00 0.00 


0.00 0.00 0.00 
0.00 0.00 0.00 
0.07 0.07 0.03 
0.01 0.01 0.00 


0.00 0.00 0.00 
0.00 0.00 0.00 
0.05 0.05 0.02 
0.00 0.01 0.00 


0.06 0.03 
0.07 0.03 
0.29 0.23 
0.23 0.11 


0.05 0.05 0.05 0.02 
0.05 0.05 0.05 0.02 
0.25 0.24 0.23 0.15 
0.14 0.14 0.12 0.06 


M24 
M70 
M160 


1.17 0.02 0.00 
NC 1.15 0.20 
NC NC 1.14 


2.85 0.75 0.16 

NC NC NC 
NC NC NC 


0.06 0.04 0.01 
2.14 0.66 0.34 
NC NC 2.81 


NC 5.14 
NC NC 
NC NC 


3.28 1.83 1.11 0.32 

NC NC NC NC 
NC NC NC NC 


P70 
PlOO 


0.83 0.01 0.01 
2.16 0.01 0.01 
NC 0.11 0.02 


1.18 0.32 0.00 
4.29 1.17 0.04 
NC NC 1.16 


0.05 0.05 0.02 
0.07 0.06 0.03 
0.24 0.14 0.06 


5.74 2.68 
NC NC 
NC NC 


1.48 0.77 0.48 0.16 
4.82 2.56 1.75 0.24 
NC NC NC 1.10 


S250 
S350 
S500 


NC 0.61 0.02 
NC 6.08 0.08 
NC NC 0.47 


NC NC NC 
NC NC NC 
NC NC NC 


1.15 0.21 0.06 
NC 1.15 0.14 
NC NC 1.14 


NC NC 
NC NC 
NC NC 


NC NC NC NC 
NC NC NC NC 
NC NC NC NC 


FUV 
NUV 


0.15 0.00 0.00 
0.33 0.00 0.00 


0.10 0.01 0.00 
0.47 0.07 0.00 


0.00 0.00 0.00 
0.00 0.00 0.00 


1.18 0.30 
3.35 1.18 


0.21 0.11 0.07 0.03 
0.60 0.26 0.12 0.03 


W3.4 
W4.6 
W11.6 
W22.1 


0.60 0.00 0.00 
1.27 0.00 0.00 
2.06 0.02 0.00 
NC 0.39 0.01 


1.05 0.21 0.00 
2.19 0.56 0.00 
3.14 1.14 0.22 
NC NC 1.65 


0.01 0.00 0.00 
0.01 0.00 0.00 
0.07 0.02 0.00 
0.60 0.28 0.04 


5.02 2.37 
NC 4.29 
NC 5.70 

NC NC 


1.18 0.53 0.29 0.01 
2.43 1.18 0.83 0.04 
3.52 1.95 1.17 0.35 
NC NC NC 1.16 



Notes.- W- = is the integral of the negative values for each kernel [see eq. ([21]) ]. 
We are abbreviating I, M, P, S, GAL, and W 

for IRAC, MIPS, PACS, SPIRE, GALEX, and WISE respectively. 

NC stands for not computed: the kernel performance would be too poor. 



to be quite safe. A kernel with W- ^ 1 is somewhat aggressive in moving power around, 
but remember that self- kernels K{A =^ A} also have W- ~ 1.15. We consider kernels with 
~ 1 to be reasonable to use, with inspection of the before and after images recommended 
in regions with large gradients or bright sources. We recommend against using convolution 
kernels with W- > 1.2, as these are, in effect, attempts to deconvolve the image to higher 
resolution, with attendant risks. 

6. Gaussian PSFs 

Gaussian PSFs are commonly used in radio astronomy. It is desirable for radio tele- 
scopes to have PSFs without extended structure to avoid sidelobe contamination. The 
illumination pattern of a single dish is often designed to return an approximately Gaussian 
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50 100 150 200 

9 (arcsec) 



Fig. 6. — Performance of the kernel J^lMieO =^ S500}. With W- = 2.81, convolution 
of MIPS160/im images into SPIRE500/im resolution is risky and not recommended. The 
convolved PSF differs from the target PSF by up to 6%. See the electronic edition of the 
PASP for a color version of this figure. 

PSF0 A Gaussian PSF lacks extended wings; the fraction of the power outside radius 6 
is exp (—0^/2(7^). Because real instrumental PSFs do not fall off so rapidly, a convolution 
kernel K going from a real \E'^ to a Gaussian with similar FWHM must "move" power 
from the wings of to the core of "^b- 

For a given instrumental PSF "^a, the optimal Gaussian PSF will be such that the 
FWHMfi will be close to FWHM^, with only mild filtering by the function Ja and with W- 
not too large. 

In order to determine an optimal Gaussian FWHM for a camera A, we compute convolu- 
tion kernels from A into Gaussian PSFs with FWHM in a range of possible values. For each 
candidate FWHM, we evaluate W_. We provide three possible FWHM. The first FWHM 
is obtained by requiring that ~ 0.3, giving a conservative (very safe) kernel that does 
not seek to move too much energy from the wings into the main Gaussian core, at the cost 
of having a larger FWHM (i.e., lower resolution). The second FWHM has ~ 0.5, and 
we consider it to be a good (moderate) Gaussian FWHM to use. The third FWHM has 



Interferometric arrays have complicated sidelobe responses and would not resemble Gaussian PSFs. 
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Fig. 7. — Performance of the self- kernel K{M70 =^ M70}. All of the self- kernels have 
iy_ = 1.16, and reproduce the original PSF to within 3%. See the electronic edition of the 
PASP for a color version of this figure. 

~ 1.0. Because this kernel is somewhat "aggressive" in moving energy from the PSF 
wings into the Gaussian core, it should be used with care, inspecting that the convolved 
images do not have any induced artifact. Table |6] gives the FWHM for three such Gaussian 
target PSFs for the MIPS, PACS, and SPIRE cameras. 



Table 6: Gaussian FWHM Suitable for MIPS, PACS, and SPIRE Cameras 





Actual 


Aggressive 


Gaussian 


Moderate Gaussian 


Very safe Gaussian 


Camera 


FWHM 


with W. 


. ^ 1.0 


with W- ^ 


0.5 


with W- 


^ 0.3 




(") 


FWHM (" 


) W. 


FWHM (" ) 


iy_ 


FWHM (" 


) 


MIPS 24/im 


6.5 


8.0 


1.00 


11.0 


0.49 


13.0 


0.30 


MIPS 70/im 


18.7 


22.0 


1.01 


30.0 


0.51 


37.0 


0.30 


MIPS 160/im 


38.8 


46.0 


1.01 


64.0 


0.50 


76.0 


0.30 


PACS 70/im 


5.8 


6.5 


0.84 


8.0 


0.48 


10.5 


0.31 


PACS 100/im 


7.1 


7.5 


1.10 


9.0 


0.52 


12.5 


0.31 


PACS 160/im 


11.2 


12.0 


1.05 


14 


0.50 


18.0 


0.33 


SPIRE 250/im 


18.2 


19.0 


1.05 


21.0 


0.44 


22.0 


0.30 


SPIRE 350/im 


25.0 


26.0 


0.98 


28.0 


0.50 


30.0 


0.27 


SPIRE 500/im 


36.4 


38.0 


0.96 


41.0 


0.48 


43.0 


0.30 
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As an example, Figure |8] shows the performance of the kernels for the PACS 160/im camera 
going into Gaussian PSFs with FWHM = 12" (aggressive, W_ = 1.05 D = 0.05), 14" (moderate, 
W_ = 0.50, D = 0.02), and 18" (very safe, W_ = 0.33, D = 0.01). The left panels are the 
performance analyses, similar to those of Figure H] - [71 The right panels show the kernel 
Fourier transform FT{K) (equation (fT9l) . In the right panel we include the unfiltered kernel 
Fourier transform (= FT0(\Ef Gauss) /-^^^(^pieo)) for comparison (dotted line). We observe 
that the filtering is only important in the Gaussian kernels with (small) FWHM = 12" and 
14". 
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Fig. 8. — Performance of the kernels K{P160 =^ Gaussian} (left panel), and their FT 
(right panel). We show the kernel Fourier transform FT{K) (equation f|T9l) ). and also the 
unfiltered kernel Fourier transform (= FT^{'^ b) / F'T^l'^ a)) ■ The filter /pieo has little impact 
on K{F1Q0 GaussianlS" } {W_ = 0.33), has moderate effect on K{F1Q0 Gaussianl4" } 
{W- = 0.50), and large effect on K{F1Q0 Gaussianl2" }(H/_ = 1.05). See the electronic 
edition of the PASP for a color version of this figure. 
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7. Usage of the Kernels 

The kernels K{A ^ B} computed here are given on a 0.20" grid. Before performing 
an image convolution, the kernel K{A =^ B} should be resampled onto a grid with the 
same pixel size as the original image I a- The resampled kernels should be centered (to 
avoid shifts in the image) and normalized so that J J K{A =^ B}{x,y)dxdy = 1 to ensure 
flux conservation. The flux in the image to be convolved should be in surface brightness 
units. After convolving the image Ia with the kernel K{A =^ B}, the resulting image will 
be expressed in the original image grid and original surface brightness units, but with PSF 

Table |5] also summarizes the kernels available. For each camera A we construct all 
the kernels K{A =^ B} with FWHMb > FWHMa/1.35 (i.e., the kernels that degrade 
the resolution or sharpen it up to 35%) plus the self- kernels K{A =^ A}. Kernels with 
FWHM^ > FWHM^ (that have larger values) tend to perform poorly and should be 
used with care. We do not recommend using any kernel with ^1-2 

As an example of the performance of the kernels applied to real (noisy) images, in 
Figure [9] we show the result of convolving images of the barred spiral galaxy NGC 1097 (after 
subtraction of a "tilted-plane" background from each image)into the S PIRE 250 Atm P SF. The 



PACS images have been reduced using the Scanamorphos pipeline (IRoussell l201ll ). Visual 
inspection of the images in Figure [9] shows them to be very similar in morphology; the 
convolution does not appear to have introduced any noticeable artifact. 

Figure ITU] shows the results of convolving the SPIRE 250 /^m image into different PSFs. 
The top row (left) shows the original^ image, (center) the image convolved with K{S250 =^ 
S250}, and (right) the image convolved with a very aggressive kernel into a Gaussian 
PSF with FWHM=18" (iy_ = 1.47). The bottom row shows the image convolved to 
the recommended Gaussian PSFs: (left) FWHM = 19" {W^ = 1.05), (center) FWHM = 
21" {W^ = 0.44), and (right) FWHM = 22" {W_ = 0.30). As discussed in §5, we recom- 
mend against using kernels with > 1.2. Visual inspection of the upper-right image in 
fact shows artifacts where the kernel (with W- = 1.47) has moved too much power out of 
some interarm pixels, which have been brought down to unreasonably low intensities. In the 
convolutions to the suitable Gaussian PSFs (bottom row in Fig. [TOj) energy is moved from 
the interarm regions into the bright nucleus and spiral arms, but the intensity levels in the 
interarm regions seem reasonable. The power that is removed from the interarm regions was 



^^By "original" image we refer to the SPIRE 250^m image delivered by the HIPE pipeline, with subsequent 
subtraction of a "tiltcd-planc" background. This is the original image that is convolved to produce the other 
images in Fig. [101 
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presumably originally power from the arms that was transferred by the wings of the SPIRE 
250/imPSF. 



Using the kernels described in the current work, lAniano et al.l (1201 if ) studied resolved 
dust modeling for NGC 628 and NGC 6946, two galaxies in the KINGFISH galaxy sample 
( iKennicutt fc KINGFISHI 120111 ). using images obtained with Spitzer and Herschel Space 
Observatory. 



8. Summary 

We present the construction and analysis of convolution kernels, for transforming images 
into a common PSF. They allow generation of a multiwavelength image cube with a common 
PSF, preserving the colors of the regions imaged. 

We generate and make available a library of convolution kernels for the cameras of the 
Spitzer, Herschel Space Observatory, GALEX, WISE, ground-based telescopes, and Gaus- 
sian PSFs. All the kernels are constructed with the best PSF characterizations available, 
approximated by rotationally symmetric functions. Deviations of the actual PSF from cir- 
cular symmetry are characterized by an asymmetry parameter given in Table [U Table 
|5] summarizes the kernels available and their negative integral a measure of their per- 
formance. We recommend using only kernels with W- ^1-2. In Table [6] we give a set of 
optimal Gaussian FWHM that are compatible with MIPS, PACS, and SPIRE cameras. 

All the kernels, IDL routines to use the kernels and IDL routines to make new kernels, 
along with detailed analysis ot the generated kernels, are puUidy a»ilable0 

We thank Roc Cutri and Edward Wright for their help providing the Wide-field Infrared 
Survey Explorer point-spread functions (PSFs); Markus Nielbock for his help providing the 
Photocamera Array Camera and Spectrometer for Herschel PSFs; and Richard Bamler, 
James Gunn, Robert Lupton and the anonymous referee for helpful suggestions. 

See the electronic edition of the PASP for a color version of this figure. 

This research was supported in part by NSF grant AST- 1008570 and JPL grant 1373687. 



^^See http://www.astro.princeton.edu/~draine/Kernels.html Kernels for additional cameras and updates 
will be included when new PSF characterizations become available. 
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Fig. 9. — Spitzer and Herschel images of NGC 1097 convolved to a SPIRE 250yumPSF. The 
SPIRE 250/im camera was convolved with the kernel i^{S250 =^ S250}. The color bar has 
the same dynamic range (lO^'^) for all images. See the electronic edition of the PASP for a 
color version of this figure. 
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Fig. 10.— SPIRE 250/im image of NGC 1097. Top row left: Original SPIRE image. Center: 
Image convolved with i^{S250 =^ S250}. Right: image convolved into an extremely aggres- 
sive Gaussian PSF with FWHM = 18" [W^ = 1.47). Bottom row: Image convolved into 
suitable Gaussian PSFs. Left: (aggressive) FWHM = 19" (W^ = 1.05). Center: (moderate) 
FWHM = 21" {W- = 0.44). Right: (very safe) FWHM = 22" {W- = 0.30). All the images 
have the same color bar. See the electronic edition of the PASP for a color version of this 
figure. 
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